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ABSTRACT 

We report on an observation of SGR 1627—41 made with the Chandra X- 
ray Observatory on 2011 June 16. Approximately three years after its outburst 
activity in 2008, the source's flux has been declining, as it approaches its quiescent 
state. For an assumed power-law spectrum, we find that the absorbed 2-10 
keV flux for the source is 1.0i J x 10~ 13 erg cm -2 s -1 with a photon index 
of 2.9 ± 0.8 (Nh = 1.0 x 10 23 cm" 2 ). This flux is approximately consistent 
with that measured at the same time after the source's outburst in 1998. With 
measurements spanning 3 years after the 2008 outburst, we analyze the long-term 
flux and spectral evolution of the source. The flux evolution is well described 
by a double exponential with decay times of 0.5 ± 0.1 and 59 ± 6 days, and 
a thermal cooling model fit suggests that SGR 1627—41 may have a hot core 
(T c ~ 2 x 10 8 K). We find no clear correlation between flux and spectral hardness 
as found in other magnetars. We consider the quiescent X-ray luminosities of 
magnetars and the subset of rotation-powered pulsars with high magnetic fields 
{B > 10 13 G) in relation to their spin-inferred surface magnetic-field strength, 
and find a possible trend between the two quantities. 

Subject headings: pulsars: individual (SGR 1627—41) — stars: magnetars — 
stars: neutron — X-rays: bursts 



department of Physics, Rutherford Physics Building, McGill University, 3600 University Street, Mon- 
treal, Quebec, H3A 2T8, Canada 

2 Lorne Trottier Chair; Canada Research Chair 

3 Space Sciences Laboratory, University of California, Berkeley, CA 94720, USA 

4 Columbia Astrophysics Laboratory, Columbia University, 550 West 120th Street, New York, NY 10027, 
USA 

5 Department of Astronomy & Harvard-Smithsonian Center for Astrophysics, Harvard University, 60 
Garden Street, Cambridge, MA 02138, USA 



-2- 



Introduction 



Magnetars are neutron stars with ultra-strong mag netic fields, the decay of which is 



theorized to power the observed radiation from the star ( [Thompson fc Duncan! Il995l . 11996 
Thompson et al.l 120021 ). Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars 



(AXPs) are two observational manifestations of magnetars. The former show repeated soft 
gamma-ray bursting activity and have relatively hard spectra, and the latter have soft spectra 
typically characterized by a blackbod y plus power law. However, the distinction between the 



two has become increasing ly blurred ( IThompson et al.l 120021 ; iGavriil et al.l |2002| ; iPons et al. 



20071 ; iPerna fc Pond 120 111 ). Magnetars generically show X-ray outbursts, which are sudden 
increases of luminosity by orders of magnitude for days to months. During outbursts, al- 
most ^jdl_J^he_gro^^ change 



see 



Woods fc Thompson! 12009 : iKaspil 120071 : Mereghettil 120081 : iRea fc Espositoll201ll . for re- 



views). The X-ray luminosities of the first-discovered magnetars were typically ~ 10 erg 
s _1 , orders of magnitude greater than their spin-down luminosity. However, more recent 
discoveries of magnetars in outburst suggest most magnetars in quiesc ence may be far less 



lumin ous (e.g. IE 1547.0-5408, XTE J1810-19 7, Swift J1822-1606; iGelfand & Gaensler 



2007 ; Bernardini et al. 2009 ; Scholz et al. 2012 ). The spin periods of magnet ara3 



are m 

the narrow range of 2-12 s and are relatively long compared to those of radio pulsars. The 
magnetic-field strength inferred from the spin period and spin-down rate is typically B > 10 14 

e formula for magnetic braking, although several wit h 



G, assuming the standard vacuum dipo 
B < 10 14 G have recently been found (IRea et al.l 12010 



Scholz et al. 



20121 ; IRea et al.ll2012h . 



The recent discoveries of a magnetar- like outb urst from a high-5 rotation-powered pul- 
sar (RPP) (PSR J 1846— 0258, Gavriil et al.ll2008l ). pulsed radio emission from a magnetar 
(XTE J1810-1 97, ICamilo et al-lbood ). a low magnetic field magnetar (SGR 0418-5729, 
Rea et al.llioioh . and the very low X-ray luminosities measured for several magnetars are of 
particular interest. These raise important questions about the relationship between magne- 
tars and other types of neutron stars, and why there is such an apparent diversity in magnetar 
properties. For example, what determin es a magnetar's quiescent X- r ay luminosity? Accord 



ingto conventional magnetar models (IThompson fc Duncan! Il995l . Il996l ; IThompson et al. 



20021 ) . there should be a correlation between B and the X-ray luminosity. However, transient 
magnetars, with their faint quiescent luminosities compared with non-transient magnetars of 
the same inferred B, challenge this. And relatedly, are high- 1? RPPs in general m agnetars in 



quiescence? PSR J1846— 0258's outburst suggests this, but it is only one source. IPons et al. 



1 See the online magnetar catalog for 



compilation of known magnetar properties, 



http://www.physics.mcgill.ca/^pulsar/magnetar/main.html| 
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( 120071 ) and lPerna fc Pond (120111 ) argue that magnetars and RPPs are all related on the basis 
of magneto-thermal evolution theory. Also a possible connection between high- 1? RPPs and 
magnetars has been sugg ested on the basis of possib l y high thermal temperatures of high- 
B RPP X-ray emission faaspi & McLaughlin] bood : Izhu et al.l bood : blausen et aD boiol ; 
Zhu et al.l 120 lit iNg fc Kaspill201ll ). However, it is important to study many magnetars and 
RPPs in quiescence to better address this question. 

SGR 1627—41 was discovered on 1998 June 15 with the Burst and Transient Source 



Exper iment (BATSE, lFishmanlll989l ). It was then identified as a SGR by iKouveliotou et al. 
dl998h . It is located at R.A. = 16 h 35 m 51 s .844, Dec. = -47°35'23".31 (J2000.0) and is 
estimated to be 11.0 ± 0.3 kpc away based on an apparent association with a star- fo rming 
region and molecular cloud (IHurley et al.lll999l ; ICorbel et al.lll999t IWachter et al.ll2004[ ) . The 
spin period and the spin-down rate were not measured until recently due to the faint nature 
of the source in quiescence. After another outburst in 2008 May, the spin period and the 
spin-down rate were measured to be 2.594578(6) s and 1.9(4) x 10~ n s s _1 , which imply an 
i nferred surface dipola r magnetic-field strength of B = 3.2 x 10 19 (PPp 2 G = 2 x 10 14 G 
jEsposito et al.ll2009al lbh. The lowest flux ever measured f or this magnetar is 6 x 10 erg 



cm -2 s~ x (2-10 keV, ~10 years after the 1998 outburst, lEsposito et al.l 120081 ) . However, 
whether it was in quiescence at that time is not clear; the luminosity might have declined 
further had the outburst not occurred in 2008. 

Here, we report the results of an observation made with Chandra in the direction of 
SGR 1627—41. We then combine our flux measurement of SGR 1627—41 with previous 
values to determine the long-term flux evolutions after its outbursts, and attempt to fit 
them to a model to infer thermal properties of the source. We compare the flux we measure 
with that from the same time after the first outburst and the lowest flux ever measured. 
Finally, we compile quiescent X-ray luminosities of magnetars and high-5 RPPs to search 
for a correlation with spin-inferred B, as might be expected in the magnetar model. 



2. Observations 



SGR 1627—41 was observed with Chandra as part of a Norma arm survey (PI: Tomsick) 
on 2011 June 16. Data from SGR 1627—41 were recorded in two consecutive 19 ks expo- 
sures (IDs 12528, 12 529), using the Advanced CCD Imaging Spectrometer I-array (ACIS-I, 
Garmire et al.ll2008l ). The data were initially processed at the Chandra X-ray Center (CXC) 
with ASCDS Version 8.3.4. After obtaining the data from the CXC, we performed all subse- 
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quent processing with the CIAO 4.3.1 softward^. We used the CIAO program chandra_repro 
along with CALDB 4.4.5 to produce the "level 2" event lists that we used for further analysis. 



3. Data analysis and results 

3.1. Imaging Analysis 

We detect the source in both exposures with wavdetect in CIAO. The best positions 
obtained are R.A. = 16 h 35 m 51 s .882, Dec. = -47°35'23".35 (J2000.0) with off-axis angle of 
7M5 for exposure 1 (ID 12528) and R.A. = 16 h 35 m 51 s .786, Dec. = -47°35'22".35 (J2000.0), 
with off-axis angle of 7'. 60 for exposure 2 (ID 12529). On average, the position is R.A. 
= 16 h 35 m 5 1 s 834 and Dec. = -47°35'22".85 (J2000.0). Using the empirical formula of 



Hong et al.l (120051 ). we estimated the position uncertainties to be P err = I'M (ID 12528) 
and P err = 1".3 (ID 12529). Therefore, we conclude the positions measured from the two 
exposures are consistent with each other and with the known position of the source. 

After finding the source location, we conducted an imaging simulation using Chandra 
Ray Tracei§ (ChaRT) and the MARJ0 tool in CIAO 4.3., and comparing with our data, find 
that there is no extended emission with 90% confidence. We find no evidence for the diffuse 



emission and hard-spectrum 'spot' reported by lEsposito et al.l (j2009al ) on the basis of a deep 



XMM-Newton observation, although this is not unexpected due to the low brightness of the 
diffuse emission (~ 10" 14 erg cm" 2 s _1 arcmin -2 ) and the relatively large angular separation 
of the hard-spectrum 'spot' (~2' from the source). 



3.2. Timing Analysis 

Since the timing resolution of ACIS-I is 3.24 s, it is not possible to measure the spin 
period of the source (~2.6 s) with this observation. However, to search for any type of 
aperiodic variability in the emission, we produced light curves with different binnings (2-7 
bins over 38 ks, to have the average counts per bin greater than 20). We also performed the 
Kolmogorov-Smirnov test (KS-test) and the Gregory- Loredo test (using the glvary tool in 
CIAO) on the data with a null hypothesis that the events are drawn from a flat light curve 



2 http://cxc. harvard . edu/ciao-4.3 
3 http:/ /cxc. harvard.edu/chart/ 
4 http: / /space.mit.edu/CXC /MARX 
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and find that the null hypothesis is statistically consistent with the data. Therefore, we 
conclude that there was no aperiodic variability in this observation over the range searched. 



3.3. Spectral Analysis 

We extracted the source counts using a radius of 10" (which includes ~95% energy 
with the spectrum given below), and the background with an annulus of inner radius 10" 
and outer radius 60" centered at the source positions reported in Section 13.11 and produced 
spectra using the SPECEXTRACT tool of CIAO 4.3 with CALDB 4.4.6.1. We find 65 ± 8 and 
60 ± 8 source counts in 2-8 keV for exposures 1 and 2, respectively. The 2-8 keV band 
is chosen because there is only 1 event below 2 keV in the source region for each exposure 
due to the large hydrogen column density (N H ). Since there are very few counts in each 
exposure, we combine the spectra for the two. Even summing the spectra, we are not able 
to tell whether a power law or a blackbody fits the data better. Therefore, we report the 
results from fits with both models. We use the wabs*powerlaw and the wabs*bbody models 
of XSPEC 12.7.C0 to fit the data. 

We fit the data with Churazov weighting with XSPEC because of the low count rate. In 
fitting, we fix Nh to the previously mea sured value, determined when the source was bright 



(1.0 x 10 23 cm 2 , lEsposito et al.ll2009af ). From the fit, we find the power-law index (r) of 
2.9 ± 0.8 (kT = 0.85±g;?| keV), and observed flux F x = l.OtgJ x lO" 13 (0.75t° f 7 x 10~ 13 
erg cm -2 s _1 ) in the 2-10 keV band for the power law (blackbody). We present the fit 
results in Table [U We also tried alternative methods to fit the data: the usual Chi-squared 
fit (grouping 15 counts per bin), a C-statistic fit (CSTAT in XSPEC, unbinned), and Chi- 
squared fit with Gehrels weighting (grouping 5 counts per bin), and obtained consistent 
results. 



3.4. Flux and Spectral Index Evolution 



Th e flux evolution of t his magnetar after the t wo outbursts was considered in several 



studi es (TWoods et al. 



1999 



Kouveliotou et al.ll2003t Mereghetti et al.ll2006at lEsposito et al. 



20081 ). iKouveliotou et al.l ( 120031 ) tried to explain an apparent plateau in the observed flux 
evolution of the 1998 o utburst (bet ween days 400 and 800) using a crustal cooling model 
( ILyubarsky et al.l 120021 ). However, iMereghetti et al.l (j2006al ) argued that the unabsorbed 



http: / /heasarc. gsfc.nasa.gov / docs / xandau/xspec 
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flux does not show a plateau and claimed that the dat a out to 1000 days w ere fitted by a 
power-law decay with index of 0.6 (F(t) oc (t — to)" ' 6 )- lEsposito et al.l (120081 ) explained the 
flux decay of the 2008 outburst (for ~1 month from the outburst) with a steep phase and a 
shallow phase. 

We show the flux ev olutio n curves includ i ng ou r measu rements and previous ones from 
Mereghetti et al.l ( ]2006ai ) and lEsposito et al.l (120081 l2009al Jbl) and show them in Figure [TJ 
We tried to fit the data out to ~3000 days with a power law (F(t) = Fi(t — t )~ a ; adding a 
constant quiescent flux did not improve the fit) or an exponential (F(t) = F\e~^ t ~ t °'' Ts - + Fq) 
or a double exponential (F(t) = Fie~ ( *"* o)/ri + F 2 e- ( *-* o)/r2 + Fq), but did not obtain a 
satisfactory fit (reduced x 2 °f 8-3, 3.0 and 3.9, respectively) for the 1998 outburst. Also a 
single power law does not describe the flux evolution after the 1998 outburst when we include 
all the data out to 10 years. 

We did not obtain a good fit for the flux evolution of the 2008 outburst with a single 
component fit. However, the data are well fit by a double exponential with decay constants 
of 0.5 ± 0.1 and 59 ± 6 days (reduced x 2= 0-7). The fit results are summarized in Table HJ 
This relaxation trend after the 2008 outburst is similar to what has been observed for other 



magn et ars (e.g. IWoods et al.l l2004j ; iRea et al.l l2009t iGavriil et al.l l2008t iLivingstone et al. 



201ll ). Note that the individual Swift XRT observations (which were in PC mode) after the 
2008 outburst did not have enough counts for a meaningful spectral analysis, and thus all 
the observat ions were assumed t o have identical spectra in determining the fluxes plotted 
in Figure [TJ (lEsposito et al.l 120081 ) . If th e spectral index changed significantly over the first 
30 days as seen in other magnetars (e.g. IZhu et al.ll2008l ). the flux values in Figure [TJ would 
have changed too, but only marginally for reasonable assumption^. 

We also plot the spectral index evolution curves for the same data in Figure [21 in order 



(Rea et al. 


2005 


Woods et al. 


2007; 


CaniDana et al. 


2007; 


Tarn et al. 


2008 


Scholz & Kaspi 


2011 


). It seems that there is no such correlation in this 



1998 outburst and only marginal evidence for it in the 2008 outburst. Also, the source 
spectrum was significantly softer just after the outburst in 1998 than in 2008. To see if 
there is any bias in spectral parameters caused by the large point spread function (PSF) of 
BeppoSAX and ASCA (> 2'), using 4 Chandra observations (ID: 1981, 3877, 12528, 12529), 
we looked for any strong soft source that could have affected the low-energy flux in the 



6 If we assume the same T vs Lx relation as Ar 
expect the first Swift data point to go down by ~10% and the rest to go up by ~10%, making t\ 
and T2 = 58 ± 6 days. 



O.lAFz (1CT 11 erg cm" 2 s~M (|Zhu et al.ll2008l ). we 

0.6±0.1 
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BeppoSAX and ^IS'CA observations. However, we find no appropriate source in a radius 
of 2'. We also checke d cross calibration results of different X-ray satellites done with Crab 
Nebula observations ( IKirsch et al.ll2005l ) but the difference in spectral indices as reported by 
different instruments for SGR 1627—41 appears too large to be due to instrument calibration 
issues alone. 



4. Discussion 

We have measured the spectrum and flux of SGR 1627—41 three years after its 2008 
outburst and find them to be consistent with those measured at the same interval following 
the 1998 outburst. We also compile the fluxes measured at different times after the two 
outbursts. The cooling curve after the 2008 outburst is well fitted by a double exponential, 
while we do not obtain an acceptable fit to the cooling curve after the 1998 outburst. We 
find at most marginal evidence for a hardness/flux correlation following either outburst. 



4.1. Flux and Spectral Index Evolution 



The flux we measured for SGR 1627—41 in 2011 June is marginally higher (by 1.5 er) 
compared to the lowest flux previously measured for this magnetar, suggesting that at this 
flux, the source is near or at quiescence. However, if we compare the spectrum in 2011 to those 
of other magnetars having low quiescent luminosities , we find that SGR 1627— 41 is signifi- 
cantl y harder (e.g. XTE J1810— 197, IE 1547.0- 5408 . iGotthelf etHlboO^lGelfand fc Gae^sler 
20071 ). This may be because SGR 1627—41 has not yet reached quiescence. It is also possi- 
ble that the sour ce has a significant hard X-ray component above ~ 10 keV as observed in 
other magnetars ( IKuiper et al.ll2006l ); such a hard X-ray component could bias the soft X-ray 
spectrum. Indeed, the source spectrum becomes softer (larger photon index) as we lower 
the high energy bound of the spectral fit, although the result is not statistically significant 
due to large uncertainties. It will be interesting to continue to monitor the source's flux and 
spectrum to see if the flux drops even lower and to see if the spectrum beco mes softer. Also 



obser ving; the source with future hard X-ray observatories such as NuSTAR ( Harrison et al. 
20101 ) can help us to determine the hard spectral component above ~ 10 keV, although it 
may be difficult if the source becomes fainter. 

Comparing the flux and spectral evolution after the outburst, we note that at ~100 
days after the 2008 outburst, the flux was an order of magnitude lower than at the same 
time after the 1998 outburst, although they became similar after ~1000 days. Thus, the 
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functional forms of the flux decays differ, as do the spectral evolutions (see Figs. [TJand[2]). 
This suggests that some aspects of the mechanism of the 2008 outbursts and/or the flux 
decay may have been different from that of the 1998 outburst. 



Thompson et al.l ( 120021 ) suggest two different mechanisms for magnetar outbursts. One 
is a sudden change in the internal magnetic field, causing a fracture of the crust, and the 
other is a sudden relaxation of the external fields. The X-ray spectrum gets harder imme- 
diately after the burst activity in the former case as the fracture will shear magnetic fields, 
generating more currents and thus more resonant up-scattering in the magnetosphere, while 
the spectrum gets softer in the latter case as the twist of external m agnetic fields is relaxed. 
M ore detailed studi es have been conducted by lPerna fc Pond (120111 ) for the crustal effect and 
by iLyutikov I (120031 ) for the magnetospheric effect, where the former studied the evolution of 
magnetic stresses in the crust to calculate properties of outbursts such as energy distribution, 
outburst waiting time and location of starquakes, and the latter estimated the time scale of 
explosive magnetospheric reconnection events that can cause magnetar outbursts. 

The lack of measurements for the pre-burst spectrum of SGR 1627—41 makes it difficult 
to diagnose the nature of the 1998 outburst in the magnetar model, as it is not known 
whether the spectrum was softer or harder before the outburst. However, the pre-burst 
spectrum was measured for the 2008 outburst and was significantly softer than immediately 
after the outburst. Thus, the 2008 outburst was likely initiated by a crustal fracture which 
twisted the external magnetic fields. 

Flux evolutions after outbursts may be explained by the untwisting magnetic field model 



(Bel oborodov fe Thompson 



2007 



Be 



jorodovl 120091 ) and/or the crustal cooling model (e.g. 



Lyubarsky et al.ll2002l ). iBeloborodovl (120091 ) explains that sudden crustal motion can twist 
the magnetic field and eject currents into the magnetosphere. The currents are then gradually 
drawn into the star and the magnetic fields are untwisted. The energy is dissipated by Ohmic 
processes in current-carrying field lines and a large fraction of the power may be radiated 
at the footpoints of the current. This process can be strongly non- uniform and produce 
complicated flux evolutions. An exa mple of a cooling cury e in the case of a localized starquake 
(ring- twist) is shown in Figure 10 of IBeloborodovl ( 120091 ) . 



A different model by iLvubarskv et al.l (120021) explains the evolution as an afterglow 
of the crustal heating. iKouveliotou et al.l (120031 ) explained the observed flux evolution of 
SGR 1627—41 after its 1998 outburst based on this model with the assumption of a heated 
inner crust and cool core, which is characterized by a t hree-phase flux de c ay: a f ast decay 
followed by a plateau and another rapid decay. However, iMereghetti et al.l (j2006al ) reported 
that the absorption-corrected flux evolution after the 1998 outburst did not have the plateau 
which is a characteristic of inner crustal cooling (due to the large heat capacity of the inner 
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crust), and we agree. Although we do not observe a plateau in the 1998 data, we note that 
ther e might be one missed in the 1998 data before ~50 days, which is different from the one 
that iKouveliotou et al.l ( 120031 ) claimed. In this case, their model may fit the data, but with 
a different set of parameters. On the other hand, we do seem to observe a three-phase decay 
following the 2008 event. 

In Figure [3j we compare the 2-10 keV luminosity decays following the 1998 and 2008 
outbursts with models of crust cooling. To calculate the cooling curves, we follow the thermal 
evolution of the c rust after a rapid deposition of energy at the start of the outburst (see 
Scholz et al.l 12012k for a recent application to Swift J1822— 1606). W e do this by solving 
the th ermal diffusion equation with a method similar to that used by iBrown fc Cumming 
( 120091 ) to model cooling of transiently accreting neutron stars, but with updated microphysics 
(Cumming et al. 2012, in pr eparation) to acc ount for the effect of the strong magnetic field 
on the th ermal conductivity (IPotekhin 1119991) and using a T e ff-T int relation for a magnetized 
envelope (IPotekhin fc Yakovlev Il200ll ). We include both phonon and impurity scattering in 
the electron thermal conductivity, and in particular set the impurity parameter Qi mp = 3. 
We take the neutron star mass and radius to be M = 1.3 M and R = 12 km, and take the 
magnetic field strength to be 2 x 10 14 G as inferred from the spin down. 

Our calculations are in ID, but we take into account the effect of the magnetic field on 
the transport of heat by assuming a dipole geometry for the magnetic field and taking an 
average over spherical shel ls (following Potekhin & Yakovlev 2001, based on the approach of 
Greenstein &: Hartkdll983l ). This means that in our calculation, we assume that the heating 
occurs in a shell over the entire surface of the star. In reality, the heating is likely localized on 
the stellar surface, which would reduce the overall luminosity because of the smaller emitting 
area. This should not change the shape of the cooling curve significantly since the thermal 
time to the surface in the thin crust is much shorter than the timescale for lateral transport 
of heat. 



We attempt to match the observed cooling curves by choosing the amount of energy de- 
posited in the crust, its location in the crust, and the neutron star core temperature T c . The 
dependence of the energy deposition as a function of d epth is not known for m agnetar out- 
bursts, and so we adopt the simple approach (following iLyubarsky et al.ll2002l ) of depositing 



a constant energy density E^IO 25 erg cm -3 in the crust. For B = 2 x 10 14 G, this represents 
a fraction 0.6i?25% of the magnetic energy density. For the 2008 outburst (bottom panel of 
Fig. [3]), we find that heating the crust in the density range 2 x 10 9 -3 x 10 10 g cm~ 3 with 
E 2 5 = 1.4 matches the shape of the cooling curve well. For the 1998 outburst (top panel of 
Fig. [3]), the energy required is more than ten times larger, £25 = 16, and must be deposited 
deeper in the crust, at densities 1 x 10 10 -2 x 10 11 g cm -3 to match the longer decay time 
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scale. 



For the 1998 outburst, while we can match the general shape of the cooling curve, we 
cannot reproduce both the rapid drop at t ~ 1000 days and the subsequent leveling off of 
the decay at t > 1000 days with our crust models. This agrees with the conclusions of 
Kouveliotou et al.l (120031 ) . who proposed that the rapid drop was a consequence of a cold 
core in SGR 1627—41. They found that for a core temperature T c ~ 2 x 10 7 K, for example 
as would be expected if direct URCA neutrino emission operated in the core, the inner 
crust would cool rapidly by conduction of heat into the core, leading to the observed rapid 
drop in luminosity at t w 1000 days. However, the luminosity we measure in this paper at 
more than 1000 days following the 2008 outburst requires a neutron star core temperature 
of T c pd 2 x 10 s K if it is due to thermal emission from the neutron star surface, and so 
the 2008 outburst is not consistent with a cold core. For example, we show models with 
T c = 3 x 10 7 K in Figure [3j illustrating that the cooling occurs much too rapidly to explain 
the observed luminosity. The observed flux at 1000 days is much greater than expected for 
the cold core case. 

Our model predicts that the current source flux will not decline by more than an order 
of magnitude in a time scale of years. However, we note that the source may not be in the 
quiescent state yet; the flux may decline slightly to the level of the last data point after 
the 1998 outburst. In this case, our model will require a slightly lower core temperature 
(~ 1.5 x 10 8 K). 



A hardness/flux correlation is expected in magnet ar models ([Thompson et al.l 12002 



Lvutikov fc Gavriil 



Gavriil et al. 



2008 



2006 



Rea et al. 



2009; 



Qzel fc Guverll2007r) and has bee n seen in other magnetars (e.g. 



Livingstone et al.l 120111 ) . It does not, however, seem to 



exist clearly in the two outbursts of SGR 1627—41, as seen in Figure [2J This is somewhat 
unexpected in the magnetospheric untwisting model; an increase in the plasma density and 
speed after an outburst is expected following a twist of magnetic fields due to a fracture, 
which decreases as the magnetar relaxes by untwisting the fields. The increase in plasma 
increases the up-scattering probability, hence resulting in a harder spectrum, with accompa- 
nying higher flux due to heat release from the interior event. We note that SGR 1627—41 is 
not the only magnetar that does not show the correlation. For example, in SGR 1900+14, 
after its 1998 giant flare, the flux decreased by a factor of three in ~18 months while neither 



the spectral index nor the blackbody temperature changed significantly (jTiengo et al.ll2007l ). 



In our crustal cooling model, as the flux declines, a decrease of temperature is nominally 
expected so the absence of a hardness/flux correlation is also puzzling. However, the degree 
of discrepancy is unclear as presently our model is not capable of predicting spectral hardness 
evolution quantitatively. Also, we note that changes in the spectral hardness of the source 
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(essentially kT for a crustal cooling event) may not be well represented by the power-law 
photon index. A better measure might be the blackbody temperature or the soft to hard 
band flux ratio. However, it was difficult to measure the spectrum unambiguously (especially 
in the soft band) for this magnetar due to its high absorption and its low count rate. 



4.2. Correlation between magnetic field and quiescent luminosity 



The quiescen t lumin osity of this magnetar is the lowest among SGRs as reported by 
Mereghetti et al.l (l2006al ). However, there are only four SGRs whose distance, hence lu- 
minosity, is known approximately. Given that SGRs and AXPs have similar natures, it is 
interesting to compare the luminosity of SGR 1627—41 with those of other magnetars for 
which the quantity is relatively well determined and to search for a correlation with, e.g., in- 
ferred surface dipola r magnetic-field strength. A correlation between B and Lx of m agnetars 
might be expected ([Thompson et al.l 120021 ; iPons et al.l 120071 : iPerna &: Pons! 120 111 ) because, 
at least for sources of comparable age, a higher B implies greater internal heating as well as 
stronger field twisting in the magnetosphere. 

To investigate a possible correlation between the spin-inferred surface magnetic field and 
quiescent luminosity in the X-ray band {Lx-, 2-10 keV) of magnetars, we select magnetars 
with a reasonable distance estimate from the McGill online magnetar catalog. We take the 
flux and distance values from references in the catalog. However, for sources that have a two- 
component spectrum, it was not possible to obtain 2-10 keV flux unless flux normalizations 
are given in the references. In these cases, we re-analyzed the archival data to obtain the 
source flux in the 2-10 keV band. To ensure that the luminosity is in quiescence, we verified 
that the measurement was done long before/after the activity of the magnetar. The data 
are shown in Table [3] and plotted in Figure HI 

A possible trend between B and Lx of magnetars is apparent in Figure HJ To see if 
the trend is significant, we calculate Pearson's correlation coefficient (r, calculated in log- 
log scale) and Spearman's rank order correlation coefficient (r s ). With the magnetars only 
(including two candidates), we obtain r = 0.63 (r s = 0.72) with a sample size of N = 16, 
corresponding to a null-hypothesis probability of p ~ 0.004 (p ~ 0.001 for r s , 1-sided). These 
values suggest a real correlation. However, in neither test are the uncertainties on the values 
taken into account. 

In order to see the effect of uncertainties in the magnetic fields and the luminosities, 
we performed simulations. We assumed that the uncertainties are 50% for the magnetic- 
field strength, the flux and the distance, and further assumed a uniform distribution for 
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the uncertainties. With 10000 simulations, we counted the occurrences in which the null 
hypothesis could not be rejected (p > 0.05). This occurred 200/10000 (480/10000 for r s ) 
times with the magnetars (with the two candidates). Although we cannot formally reject 
the null hypothesis, this at least suggests a trend between B and the quiescent luminosity. 

For sources that have significant flux below ~ 2 keV, considering the 2-10 keV band only 
may not be optimal. Therefore, we have repeated this analysis for luminosities in the 1-4 
keV range, where thermal emission dominates but where the effects of interstellar absorption 
are more pronounced, and found similar results to those in the 2-10 keV band. 

Note that we have not included upper limit measurements in Table |3] and Figure H] 
because of the difficulty of handling them statistically in the correlation calculation and 
translating the limit to the 2-10 keV band. However, we have verified that in no case is a 
reported upper limit in clear contradiction with the observed possible trend. An important 
source to include in the future is SGR 0418+5729, as it has very low r eported field with no 
quiescent X-ray luminosity measured yet, although iTurolla et al.l (120111 ) suggest it may have 
higher-order multipoles, and thus may not lie on the trend at the spin- inferred B. 



On FigurelU we plot Lx oc B 4A , similar to the relation predicted by lThompson fc Duncan 



( 119961 ) from internal heating due to magnetic dissipation. This relation is broadly consistent 
with the data, albeit considerable scatter exists. Note that the lack of magnetars with Lx 
greater than 10 36 erg s -1 is consistent with the internal heating model whe re the X-ray lu- 
minosity saturates at Lx = 10 35 -10 36 er g s _1 due to rapid neutrino cooling ( Ivan Riper et al. 



1 99 ll : iThompson fc Duncanlll993l . Il996l ). This model, however, explains the quiescent X-ray 
flux as being thermal in origin, and thus naively would predict a correlation between the 
quiescent surface temperature and the magnetic field while there is no clear observational 
corre lation between them in the magnetar population ( jZhu et al.l l2009t iKaspi fc Boydstun 
2010|). This suggests that there can be a significant "twiste d magnetospheric" effe ct in the 
soft X-ray emission. The "twisted magnetospheric" model (IThompson et al.ll2002l ) predicts 
a correlation between B and Lx, where no sim ple relation is giv e n due to the difficulty 
in estimating the "twist" (B^/Bg). Nevertheless, IThompson et al.l ( 120021 ) explain that the 
initial output of a magnetar is provided by surface heating (L x oc B AA ) and is increased by 
a modest factor due to multiple scattering. Therefore, Lx should be a strong function of B 
in this model. 



Alternatively, iPons et al.l ( 120071 . 120091 ) showed an interesting trend between the effective 
surface temperature and the magnetic field with 27 neutron stars, including both RPPs and 
magnetars, over magnetic-field range 10 12 -10 15 G. They explain the trend with the decay 
of crustal currents, where T e ff oc B 1 ! 2 is expected in a simple illustrative calculation. If 
we assume that the luminosity is from blackbody emission (i.e. that the initial output of 
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magnetars is thermal), we expect L oc B 2 in this model, which can also roughly describe the 
possible trend we find (see Fig. H]). On the other hand, pure blackbody emission is likely 
an oversimplification so more detailed modelling, such as consideration of the effects of an 
atmosphere, is warranted. Age is also an important factor in determining the luminosity in 
this model and likely for understanding the scatter in Figure HJ 



4.3. Connection to High-B RPPs 



The 2006 outburst of the young, high-5 RPP PSR J1846-0258 flGavriil et all 120081 ) 
clearly demonstrates a connection between m agnetars and high--B RPPs. Also, a m odel of 
magneto-thermal evolution in neutron stars (IPons et al.l 120071 ; iPerna fc Pons! 120111 ). moti- 
vated by the apparent correlation between the inferred magnetic field and surface tem per- 
ature over a broad range of magnetic fields, suggested a connection (see iKaspil l2010l for 
review). It is interesting to ask if a correlation between B and Lx exists in the high-5 RPP 
population, and to search for a connection to the magnetar population. 



Using lOlausen et al.l (120101 ) and other references (e.g. lZhu et aU201ll ; Kaplan fc van Kerkwijk 



201 ll ). we plot B vs Lx of high- B RPPs in Figure HJ Interestingly, we note that these appear 
roughly consistent with the possible trend noted for magnetars alone. We consider this trend 
more quantitatively. With high- 1? RPPs only, we obtain r = 0.18 (r s = 0.1), consistent with 
the null- hypothesis as one can easily see in the plot. This is not surprising, considering 
the sample size (N = 5) and uncertainties. Also note that the luminosities (2-10 keV) of 
some high-i? RPPs are hig hly uncertain, as they were measured in a lower energy band (e.g. 
Kaspi fc McLaughlin! 120051 ) and extrapolated to the 2-10 keV band. However, if we combine 
high-i? RPPs and magnetars, we obtain a better correlation of r = 0.77 (r s = 0.82) and 
p < 0.0001 (p < 0.0005 for r s , N = 21) than we do with magnetars alone. Also our simula- 
tions to investigate the effect of uncertainties (see Section I4.2[) show that the null hypothesis 
is always rejected in this case. Repeating this analysis in the 1-4 keV band yields similar 
results. Having a better correlation with high-I? RPPs and magnetars than with magnetars 
alone may suggest that high-S RPPs and magnetars share similar physical processes, which 
evolve continuously function of magnetic field. 

There is large scatter in the correlation plot. Uncertainties in estimating the true mag- 
netic field from the inferred surface dipolar magnetic field3, distance, unabsorbed flux and 
age effects are obviously possible contributors. Further, some variation is expected depending 



7 Note that some magnetars such as SGR 1900+14 and SGR 1806—20 have a large uncertainty in the 
spin-down inferred magnetic-field strength due to spin-down rate variations. 
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on the efficiency of multiple scattering of thermal photons, and radiation localization effects 
may play a role. However, AXP 4U 0142+61 and IE 2259+586 stand out as having large 
luminosities with relatively weak magnetic fields. One possible explanation is that the spin- 
down inferred magnetic field is sensitive to the dipola r component only and these magnetars 
have very strong toroidal or multipole components ([Thompson et all l2002t iPerna fc Pons 



201ll ). X-ray polarimetric observations may be able to test this idea. 



5. Conclusions 

Using Chandra observations, we have measured the spectrum and absorbed flux in the 2- 
10 keV band for SGR 1627—41 approximately 3 years after its 2008 outburst. The spectrum 
was consistent with a power law having T = 2.9 ± 0.8 (or a blackbody having kT = 0.85l ; 2 g 
keV), and the absorbed flux was l.OtgJ X 10~ 13 erg cur 2 s" 1 (0.75+o; 2 ° x 10~ 13 erg cm" 2 s" 1 
for a blackbody spectrum) in 2011 June. Although the source flux is similar to that detected 
a comparable amount of time following its 1998 outburst and is similar to the lowest yet 
seen from this source, it is unclear whether it has reached true quiescence, as its spectrum is 
significantly harder than in other magnetars in quiescence. We showed that the flux evolution 
of the source after its outburst activity in 2008 followed a double exponential with decay times 
of 0.5 + 0.1 and 59 + 6 days. Our model fitting, assuming the flux relaxation is due to crustal 
cooling, suggests that the core temperature of SGR 1627—41 is high (T c ~ 2 x 10 8 K) and 
that the energy was deposited in the outer crust (at di fferent depths) for the two outbursts. 



This is the same conclusion as for Swift J1822— 1606 (IScholz et al.ll2012l ) and may provide 
an interesting constraint on crust breaking models. We show that the 2008 activity of 
SGR 1627—41 was likely to have been initiated by a crustal fracture, causing a twist of the 
external magnetic fields. However, for this magnetar, we see no clear correlation between 
flux and spectral hardness as seen in other magnetars, which is puzzling. Finally, we find a 
possible correlation between the inferred magnetic field and the quiescent luminosity of 16 
magnetars (including two candidates). We also note that the correlation becomes stronger 
if we include high-5 RPPs, which further suggests a connection between high-5 RPPs 
and magnetars. The discovery and detailed study of more high-5 RPPs and magnetars 
in the future will help us to better understand the physical connection between these two 
populations. 
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Table 1: Summary of the spectral fit results for an absorbed power law and an absorbed 
blackbody with Churazov weighting. 



Fit Function 


N H 


T/kT 


Flux a 


x 2 /DoF 




(10 22 cm" 2 ) 


(/kcV) 


(10 -13 erg cm -2 s _1 ) 




Power Law 


10 b 


2.9 ±0.8 


1 - u — 0.2 


12.1/21 


Blackbody 


10 b 




u -'°-0.17 


14.5/21 



Fits are conducted in the 2-8 keV band. All uncertainties are at the 90% confidence level. 



a Absorbed flux in the 2-10 keV band. 

b Frozen at the value from lEsposito et al.l (2009a) 



Table 2: Fit results for the flux evolution after the 2008 outburst. 



Fit Function 


Fi a 

(10~ 12 erg cm~ 2 s _1 ) 


a/n 
(/days) 


F 2 & T 2 
(10" 12 erg cm~ 2 s _1 ) (days) 


Fq* 

(10~ 12 erg cm" 2 s _1 ) 


X 2 /DoF 


Power Law 


7.6(3) 


0.60(2) 






51/12 


Exponential 


2.9(2) 


48(5) 




0.23(4) 


308/11 


Double exponential 


28(3) 


0.5(1) 


2.1(2) 59(6) 


0.22(4) 


6/9 



a Unabsorbed flux in the 2-10 keV band. 
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Table 3: Spin-inferred surface magnetic- field strength and 2-10 keV quiescent X-ray lumi- 
nosity of magnet ars and high- 1? RPPs. a 



OvJ Li. J. V_.l_. 1 N tHIlL, 




r X 




T) i q+~ a n r'(~> 

1 ' 1 o LojII^V^ 


^X 


Ref. 




(10 14 G) 


( ~\ O fi"cr cm 


-2 -1\ 




HO 35 pre s _1> l 




C< -rj_ T-i qoo 1 fiOft 


U.oo 


U.U4 




1.0 


l,Z X 1U 


1 


IE 2259+586 


0.59 


17.7 




4.0 


3.4 x 10"" 1 


2 


CXO J164710.2-455216 


0.95 


0.14 




5 


4.2 x 10" 3 


3 


4U 0142+61 


1.3 


70.2 




3.6 


1.1 


4 


XTE J1810-197 


2.1 


0.02 




3.5 


3.5 x 10~ 4 


5 


IE 1547.0-5408 


2.2 


0.56 




3.9 


1.0 x 10~ 2 


6 


SGR 1627-41 


2.2 


0.17 




11.0 


2.5 x 10" 2 


7 


PSR J1622-4950 d 


2.8 


0.065 




9 


6.3 x 10" 3 
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See McGill SGR/AXP onlin e catalog and refer ence s therein. For pu lsar data, see | 01ausen et al.l (120101) and 
refere nce s therein. Refs: [1] JScholz et al.ll2012h [21 Jzhu et al.ll2008h f3l JlVtuno et al.l|200i [41 (bohler et al ' 
20051) [51 JGotthelf et al.l l 2004h [61 dCeffand fc Gaenslerll2007r) [71 dEsposito et al.ll2008l) [81 (iLevin et alll201 " 



[91 (rTam et al.ll2008l) [10 1 (rTiengo et al.lbooj [111 Jden Hartog et al.ll2008l) [12] tealpern fc GottheliEoi 



2007) [17] (Zhu ct al 



20101) [21] (|Zhu et al 



[131 (ITiengo et al.l 120091) [1 41 (Kumar fc San-Har bll2010l) [ 151 dMereghetti et al.ll2006bl) [16] (Esp osito et al 



201 ll) 



2009) [18] dSafi-Harb fc KumaJ20oI) [19] dMcLaughlin et al.ll2006l) [20] ddausen et al 



a Uppcr limit measurements are not included as they cannot be used in the correlation coefficient calculation. 
b The lowest unabsorbed flux ever measured for the magnetar in the 2-10 keV band. Converted with PIMMS 
or XSPEC if 2-10 keV unabsorbed flux is not given in the reference. 
c Calculated from Fx and distance. 
d Candidates. 

c Classffied as a rotating radio transient (RRAT). 
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Fig. 1. — 2-10 keV flux values and best-fit function (solid line) for the flux evolution after 
the 1998 (top) and 2008 (bottom) outbursts. For the 1998 outburst, neither a power law 
nor an exponential decay gives a satisfactory fit to the flux evolution (an exponential fit is 
shown in the plot). The 2008 outburst data are w ell fit by a double expone n tial with decay 
constants of 0.5 and 59 days. Data are taken from iMereghetti et al.l ( j2006af ). lEsposito et al. 
J2008I . l2Q09al Jbh and this work. 
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Fig. 2. — Power-law index T vs time and 2-10 keV unabsorbed flux for the 1998 and 2008 
outburst relaxations. T for two XMM-Newto n data points at ^2000 days in the left panels 
was fixed because of the low counts statistics (jMereghetti et al.ll2006al ). Due to the low count 
rate in the Swift observations, all the Swift data are combined to o btain o n ly one spectral 



i ndex in the right panels. Data are taken from iMereghetti et al.l (j2006al ). lEsposito et al. 
J2008I . l2009al lbl) and this work. 
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Fig. 3. — Comparison of crust cooling models with the 2-10 keV luminosity decays following 
the outbursts in 1998 (top panel) and 2008 (bottom panel). All of the models shown assume 
a constant energy density £"25 10 25 erg cm -3 is deposited instantaneously in the crust at 
densities p min < p < p ma x- Top panel (1998 outburst): The solid curve is for 1 x 10 10 < 
p < 2 x 10 11 g cm -3 , E25 = 16, T c = 2 x 10 8 K. The dashed curve is for a cold core, 
1 x 10 10 < p < 2 x 10 11 g cm" 3 , £25 = 18, T c = 3 x 10 7 K. The dotted curve is for energy 
deposition extending to 1 x 10 10 g cm" 3 < p, with E 2 $ = 18, T c = 3 x 10 7 K. Bottom 
panel (2008 outburst): The solid curve is for 2 x 10 9 < p < 3 x 10 10 g cm" 3 , E 25 = 1.4, 
T c = 2 x 10 8 K; the dashed curve has 2 x 10 9 < p < 3 x 10 10 g cm" 3 , E 25 = 1.7, T c = 3 x 10 7 K. 
The dotted curve is for 2 x 10 9 < p, E 25 = 1.7, T c = 3 x 10 7 K. 
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Fig. 4. — B vs Lx (2-10 keV) of magnetars and high-i? RPPs with known distance including 
two magnetar candidates (PSR J1622-4950, CXOU J17140 5 7-381031). See Tab l e B] fo r 
data. The solid line indicates the relation, Lx oc B 4A , given by lThompson fc Duncan! ( 119961 ). 
and the dashed line shows the relation, Lv oc B 2 , obtained with the kT vs B relation of 



Pons et al.l (120071 ) and an assumption of pure blackbody emission. A possible trend between 
the surface magnetic-field strength and the luminosity can be seen. Uncertainties of 50% on 
the distance and the flux are assumed. 



